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SIMILAR  SOLUTIONS  IN 
VIBRATIONAL  NONEQUILIBRIUM  NOZZLE  FLOWS 
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Wri ght-Patterson  Air  Force  Base,  Ohio 


ABSTRACT 

The  problem  of  obtaining  similar  solu¬ 
tions  in  pseudo-one-dimensional  nozzle 
flows  uith  vibrational  relaxation  is  con¬ 
sidered.  The  governing  equations  are  trans¬ 
formed  into  a  similar  form  by  using  a  new 
similarity  parameter  n  so  that  the  nonequi¬ 
librium  flow  solutions  depend  on  two  param¬ 
eters  for  a  given  kind  of  gas.  However, 
the  similar  equations  are  further  reduced 
to  a  universal  form  by  a  transformation  of 
the  independent  variable  r)  to  £  such  that 
the  similar  solutions,  for  all  combinations 
of  initial  conditions  and  for  a  defined 
family  of  nozzle  shapes,  depend  on  a  single 
parameter  X  which  are  then  presented  in  a 
single  graph.  It  is  also  shown  that  the 
equilibrium  solutions  depend  on  only  one 
parameter  £.  The  parameters  £  and  X  are 
also  the  exact,  general  correlating  param¬ 
eters.  Uith  the  present  similar  solutions 
the  flow  quantities  in  a  nozzle  arc  readily 
determined  without  resorting  to  complex 
computer  programs. 

NOMENCLATURE 

A  »  area  ratio  (a*/a;) 

a  ■  speed  of  sound 

C  ■  constant  defined  in  Eq  (8) 

0  «  constant  defined  in  Eq  (8) 

f i , f 4 , f 3  ■  functions  dzfined  in  Appendix  A 

Ho  «■  reservoir  enthalpy  (HQ  =  Mq/u^) 

h  -  static  enthalpy  (h  =  h'/Uj) 

kj.kj  =  functions  defined  in  Appendix  B 

L  -  nozzle  scale  parameter 

(L  =  r'./tan  y) 


p  -  pressure  (p  =  p'/Pj) 

R  »  molecular  gas  constant 

(R  =  R/m) 

-  nozzle  throat  radius 

SQ  >*  reservoir  entropy  (S^  =  S^/R) 

S  »  reservoir  entropy 

r  (Sr  =  S^/R)  =  15.2  for  nitrogen 

T  .  «  characteristic  dissociation 

"  temperature 

T .  ■  translational  temperature 

<rt  *  H'V 

T„  “  vibrational  temperature 

v  (tv  =  w 

*j,  »■  characteristic  dissociation 

a 

velocity  [ud  =  (RTd)*] 

U  »  velocity  (u  =  u'/Uj) 

-  velocity  defined  as  (Pq/Pq)* 

X  “  distance  along  the  nozzle  axis 

(x  =  X'/L)  \ 

e  «  vibrational  energy  (e  =  e'/RTj) 

I  ° 

e(o>)  m  equilibrium  vibrational  energy 

T  ■  vibrational  relaxation  tine 

constant 

h  ■  function  defined  as 

loge(uA/p*u*) 

Pd  ■  characteristic  dissociation 

u  density 

p  »  density  (p  =  P'/Pj) 


K  "  flow  Mach  number 

m  *■  molecular  weight 

N  «  number  of  vibrational  levels 

(N  -  Td/ev) 

(Hs),  “  defined  in  Eq  (17) 

(H$)t  “  defined  in  Eq  (2S) 

N,.  -  defined  in  Eq  (31) 

p  .  «■  characteristic  dissociation 

pressure  (pd  = 


8  •>  characteristic  vibrational 

fo  c  fi  /T  .) 

- , - 5 —  >  -  -  y '  d 

2y  “  expansion  ang  -  of  conical 

nozzle  ou  asymptotic  cone 
angle  of  hyperbolic  nozzle 

“  parameters  defined  in  Eqs  (17) 

L  (27)  and  (34)  respectively 

$  -  vibraUional  temperature  func¬ 

tion  defined  as  $  =  { 9 V/T ^ ) 

“  translational  temperature  func¬ 
tion  defined  as  $  11  ( 0 y /T^ ) 


1 


X*.X,.X 


“  parameters  defined  in  Eqs  (24). 
(28)  and  (35)  respectively 

5  -  U0  -  n) 

Subscripts 

*  “  nozzle  throat 

Q  “  reservoir 

Superscript 

i 

“  dimensional  quantity 


X.  INTRODUCTION 

The  phenomenon  of  thermal  anu/er  chem¬ 
ical  nonequilibrium  in  flowing  gases  at 
high  temperatures  is  of  fundamental  impor¬ 
tance  in  rocket  nozzles  as  well  as  in  lab¬ 
oratory  facilities  which  are  being  used  to 
simulate  the  flow  fields  over  bodies  in 
hypersonic  flight.  In  the  past  decade  a 
considerable  amount  of  theoretical  effort 
has  been  directed  toward  understanding  the 
nonequilibriun  flow  effects  in  nozzles 
under  steady  flow  conditions.  The  problem 
has  been  studied  under  the  assumption  of 
pseudo-one-dimensional,  adiabatic,  inviscid 
flow.  In  spite  of  these  simplifying  assump¬ 
tions  the  solutions  are  far  from  being  sim¬ 
ple  and  are  often  plagued  by  many  numerical 
procedural  difficulties,  a  comprehensive 
review  of  this  problem  is  presented  in 
Ref  1.  More  recently,  a  tine-dependent 
analysis2  has  been  proposed  which  circum¬ 
vents  some  of  the  numerical  difficulties 
but  retains  the  problem  of  determining  the 
flow  quantities  through  a  nunerous-s tepped 
process.  The  basic  problem  of  obtaining 
numerical  solution*  for  nozzle  flows  with 
vibrational  energy  relaxation  has  been 
studied  by  several  authors  (e.g.,  Refs  3, 

4,  5).  The  present  state-of-the-art  for 
solving  vibrationaj.  nonequilibriun  nozzle 
flows  requires  complex  computer  programs 
with  which  the  flow  variables  are  deter¬ 
mined  by  numerical  Integration  for  any 
given  initial  and  boundary  conditions.  How¬ 
ever,  this  approach  does  not  provide  suit¬ 
able  theoretical  comparisons  for  use  by  the 
experimentalist  because  of  the  many  vari- 
ab'es  Involved.  Thus,  it  is  apparent  that 
general  correlating  parameters  are  needed. 
Several  approximate  analyses6*7  have  been 
used  to  predict  the  frozen  vibrational 
temperature  In  nozzle  flows.  In  these 
analyses  approximate  correlating  param¬ 
eters  have  been  deduced  which  do  not  In¬ 
clude  ell  the  variables  of  the  problem, 
furthermore,  these  analyses  predict  only 
the  frozen  vibrational  temperature. 

From  the  preceeding  discussion  it  is 
obvious  that  suitable  similar  solutions 
to  this  problem  are  highly  desirable.  Such 
similar  solutions  would  net  only  eliminate 
the  need  for  repeated  complex  computations 
but  a  so  provide  for  the  experimentalist, 
the  badly  needed  correlating  parameters. 


In  the  present  analysis,  the  governing 
equations,  for  a  pseudo-one-dimensional 
nonequilibrium  nozzle  flow  with  vibrational 
energy  relaxation  but  no  dissociation,  are 
transformed  into  a  similar  form  by  using  a 
new  independent  variable  n .  The  definition 
of  the  similarity  parameter  r\  and  the  meth¬ 
od  of  transformation  of  the  governing  equa¬ 
tions  in  the  present  case  are  very  similar 
to  those  used  in  the  analysis  presented  by 
the  authors8  which  deals  with  the  problem 
of  di8sociational  nonequilibrium  nozzle 
flows.  It  is  shown  that  the  similar  solu¬ 
tions,  for  a  family  of  nozzle  shapes  and 
e  specified  gas,  depend  on  two  parameters, 

S0  and  ) ,  in  add*  tion  to  the  independent 
variabl*  n.  However,  the  similar  equations 
are  further  reduced  to  a  universal  form  by 
a  transformation  of  the  Independent  vari¬ 
able  r)  to  C  so  that  the  similar  solutions 
depend  on  a  single  parameter  x  with  £  as 
the  Independent  variable.  General  similar 
solutions  which  can  be  used  for  all  com¬ 
binations  of  initial  conditions  are  pre¬ 
sented  in  a  single  graph  for  nitrogen.  The 
parameters  x  and  £  arc  the  exact  and  gen¬ 
eral  correlating  parameters.  The  limiting 
colutions,  namely  frozen  and  equilibrium 
solutions,  ore  also  obtained  with  £  as  the 
independent  variable  and  these  solutions 
depend  on-  £  only.  The  upproxlmite  corre¬ 
lating  parametere  that  have  appeared  iu  the 
literature6*7  can  be  deduced  from  the 
present  general  correlating  parameter  X  and 
the  approximations  involved  in  these  anal¬ 
yses  are  discussed. 

II.  TRAKSFOF.HATIOH  OF  GOVERNING  EQUATIONS 

In  the  analysis  of  the  dissociation  non- 
equilibrium  case3  it  was  logical  to  use  the 
characteristic  dissociation  values  to  non- 
dimensionalizc  the  governing  equations 
since  the  characteristic  dissociation  den¬ 
sity  and  temperature  appear  in  the  govern¬ 
ing  equations.  The  same  values  are  also 
usea  for  nondlmensionalization  in  the 
present  case  for  the  following  reasons.  The 
governing  equations  contain  only  the  char¬ 
acteristic  vibration* I  temperature  6y  and 
there  Is  no  corresponding  characteristic 
density.  Furthermore,  it  is  shown  later 
that  tin  similar  solution*  in  the  present 
case  do  cot  depend  on  the  values  used  for 
nondloenslonsll.'.ation  and  hence  can  be 
arbitrarily  chocen. 

The  governing  aquations  for  a  steady, 
pseudo-one-dimensional,  adiabatic,  inviscid 
flow  with  negligible  dissociation  are  con¬ 
sidered.  fhe  flow  variables  d'.  o',  T’ 
and  u'  arc  nondimensionalized  with  the 
corresponding  iissociation  values  which  are 
defined  in  tnc  nomenclature.  The  nozzle 
area  and  the  distance  along  the  nozzle  axis 
are  nondimcp  ilonalized  by  the  nozzle  throat 
area  A*  and  ncale  parameter  1,  respectively. 
The  governing  equations  after  nondimen- 
slonalizotion  are  wrlttcu  as 

Equation  of  state: 

P  K  pTt  (i) 


2 


Conservation  of  mass: 


Conservation  of  energy: 


puA  = 


(2) 


Conse i va t lor.  of  energy: 


u?/2  +  h  =  Hq  =  const 


where 


(3) 


h  =  ( 7  /  2 )  T  t  +  e 

Conservation  of  momentum: 

udu  +  dp/p  =  0 

Rate  equation  (Landau-Teller  type)  : 


(4) 


de 

dx 


uudTv 


[e  (“} 


e!l 


(5) 


where  £(<")  corresponds  to  vibrational  equi¬ 
librium.  Applying  .'he  cut-off  harmonic 
oscillator  approximation  (see  Appendix  A), 
the  vibratio  energy  C  can  be  given  in 
terms  of  vibrational  temperature  as 


0  N9 

e°-/Tv-l  efi°/Tv-l 


(6) 


where  a  Boltzmann  distribution  of  the  mo¬ 
lecular  number  density  in  the  different 
vibrational  levels  is  assumed.  The  corres¬ 
ponding  vibrational  energy  at  vibrational 
equilibrium  e(“)  is  given  as 


e(«') 


0  __N0__ 

e6/Tt-l  '  eK0/Tt-l 


(7) 


The  vibrational  relaxation  time  constant  t.( 
has  been  correlated1  for  many  gas  systems 
by  using  measured  values  and  may  be  rep¬ 
resented  by  an  expression  of  the  form 


T  P 
VK 


=  Ce 


m  cv/tJ) 


i/3 


(8) 


where  C  and  1)  are  constants  for  a  given  gas. 


An  independent  variable  *1  is  now  de¬ 
fined  as  n  ’  loge (uA/p*u*)  and  the  tempera¬ 
ture  functions  are  defined  as  $  <*  0/Tv  H 
0v/Ty  and  'V  ■  0/Tt  =  8V/Tt'  Using  these 
definitions  and  £qs  (7)  and  (8),  the  gov¬ 
erning  equations  are  reduced  to 


-  +  0 


11  + 


1 


,Nf 


•1_ 


■  H. 


L  2  e*-l  o‘ 

Conservation  of  momentum: 

dn-  (5/2)(dti>/«|,}-(l-f1).|»e*d4>/(e4-'i)i 


Rate  equation: 

d$  Lopde"ne"I)'i' 


1/3 


dx  Cuud 


(11) 

0  (12) 

±!\  l±hl 

e*-l  /'  (1-f,) 

(13) 


where  f4  and  f2  are  the  factors  which  take 
into  account  the  effect  of  the  cut-off 
harmonic  oscillator  approximation  end  are 
given  in  Appendix  A.  It  is  shown  in  Appen¬ 
dix  A  that  these  f  functions  tend  to  zero 
for  »  “,  corresponding  to  the  simple 

harmonic  oscillator  approximation.  Further¬ 
more,  it  is  also  shown  in  Appendix  A  that, 
for  temperatures  even  as  high  as  7000°K  for 
nitrogen,  these  functions  are  negligible. 
Hence,  the  simple  harmonic  oscillator  ap¬ 
proximation  is  sufficiently  accurate  for 
the  present  analysis. 


The  independent  variable  x  in  Eq  (13) 
is  changed  to  n  by  using  the  following 
procedure.  From  the  conservation  of  mo¬ 
mentum  and  mass  the  following  relation  can 
be  derived:  Q4) 

d  iogep/d  1ogeA  =  -(dr,/d  logeA)  =  M2/(l-H2) 

where  M  »  u/a  and  a2  »  dp/dp.  For  a  family 
of  nozzles  with  the  area  distribution  given 
by  A  ■*  (1  +  xi)1  and  using  Eq  (14)  the  fol¬ 
lowing  expression  can  be  written: 


it.  »  /  «L\ 

dx  \Ha-l  j 


(15) 


Substituting  for  A~^7ij  from  the  definition 
of  n,  Eq  (15)  can  be  written  as 


ul/iJ  e-n/iJ  (1-A-1'1)(J“1>/J 


(16) 


Combining  Eq  (13)  and  Eq  (16)  the  rate 
equation  is  transformed  to 


Equation  ot  state: 


PV 


ee 


-n 


(?) 


d£ 

dn 


<»s>. 


2  -i 

e  e 


1/3 

-tn(  1-1/iJ  )+<•*  D+^l 


Conservation  of  mass: 


o  " 


e 


-n 


(10) 


(«*-.♦> 


(17) 
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where 


for  Y  "  7/5.  Equation  (19)  io  the  familiar 
isentropic  relation  foi  a  perfect  gas. 


X2  B  1oge  |(P*uJ1/iJ  LpdUd0/1jcJ 

and 

(»,),  ■  (h^t) 

The  problem  under  consideration  has 
been  reduced  to  solving  two  differential 
equations,  Eqo  (12)  and  (17),  for  two  un¬ 
knowns  <j>  and  <j>  with  r)  as  the  independent 
variable.  Once  p  and  $  are  determined  by 
solving  Eqs  (12)  and  (17),  the  other  un¬ 
knowns,  p,  p  and  u  are  obtained  from  the 
other  governing  equations  which  are  simple 
algebraic  equations. 


Equilibrium  Flow  Solution 

This  limiting  case  is  achieved  when  the 
vibrational  relaxation  time  tv  is  very 
short  and  this  happens  if  X2  «®.  With 
this  situation  the  condition  of  $  ■  tj)  can 
be  inferred  from  f ne  rate  equation  (Eq 
(17)).  With  <J>  ■  p,  the  generalised  momen¬ 
tum  equation  (Eq  (12))  can  be  integrated 
and  given  as 

n  -  (5/2)  loge^  -  loge( l-e-'3')  + 

^/(e^-l)  +  f3  -  const  (20) 


Character  of  Similar  Equations 

The  main  motivation  in  expressing  the 
rate  equation  in  the  form  shown  in  Eq  (17) 
is  to  combine  all  the  parameters  of  the 
problem  into  a  single  parameter  X,.  How¬ 
ever,  Eq  (17)  also  contains  the  additional 
parameters  D  and  i j .  The  parameter  D  is 
a  constant  for  a  given  gas  and  this  means 
separate  solutions  have  to  be  obtained  for 
each  kind  of  gas.  The  parameter  ij  is  the 
nozzle  shape  parameter  and  a  given  value  of 
ij  covers  a  family  of  nozzle  shapes;  its 
value  is  2.0  for  nozzle  shapes  of  practical 
interest,  namely  conical  (i  «*  2,  j  «  1)  and 
hyperbolic  (i  *•  1,  j  *»  2)  shapes.  Equation 
(17)  also  contains  a  rather  undesirable 
term  (Ns)2  which  is  a  function  of  the  flow 
quantities  in  the  nozzle  and  is  not  a  con¬ 
stant.  In  this  respect  this  term  may  be 
called  a  nonsimilar  function.  The  proper¬ 
ties  of  (Hs).  and  a  method  of  including 
the  effect  of  (Ng)2  into  the  solutions  of 
similar  Eqs  (12)  and  (17)  are  discussed  in 
one  of  the  following  sections.  In  addition 
to  specifying  X2 ,  D  and  i j ,  the  initial 
values  of  p,  for  a  given  value  of  n, 
have  to  be  specified  for  the  solutions  of 
Eqs  (12)  and  (17).  The  specification  of 
these  initial  values  will  be  facilitated  by 
the  limiting  solutions,  nai  ely  the  frozen 
and  equilibrium  solutions. 


III.  LIMITING  SOLUTION'S 

Frozen  rlow  Solution 

Tills  limiting  case  is  achieved  when  the 
vibrational  relaxation  time  Tv  is  very 
large  and  this  happens  if  X2  -  -■>.  Then 
d$/dn  •*  0  follows  from  the  rate  equation 
and  hence  the  trivial  solution  p  «  constant. 
The  generalized  momentum  equatiou  (Eq  (12)) 
simplifier  to 

dn  -  (5/2)(dv/<{-)  =  0  (18) 

This  equation  car.  be  integrated  and  given 
as 


[l/{ Y-l)  ] 


where  the  function  f}  is  a  term  associated 
with  the  cut-off  harmonic  oscillator  ap¬ 
proximation  and  is  negligible  as  shown  in 
Appendix  A.  Equation  (20)  is  the  expres¬ 
sion  for  the  change  in  entropy  in  3n  equi¬ 
librium  flow.10  This  equation  also  shows 
that  entropy  Is  conserved.  The  expression 
for  entropy  in  this  case  can  be  given  as 

So  =  (s;/R)  =  n  -  (5/2)  loge«/  - 

logeO-e'*)  +  +  U  +  Sr  (2i) 

where  Sr  «  S£/R  and  Is  the  reference  en¬ 
tropy.  The  equilibrium  flow  solution  can 
be  obtained  from  Eq  (21)  with  U  as  the  In¬ 
dependent  variable  and  SQ  as  the  parameter. 
It  is  worth  noting  that  the  limiting  solu¬ 
tions  do  not  depend  on  the  nozzle  geometry. 


NONEQUILIBRIUM  SOLUTIONS 


Initial  Values 


To  obtain  .leneqai  librium  solutions,  a 
set  of  initial  flow  quantities  if>,  <J>  and  n 
have  to  be  specified  before  solving  the 
similar  governing  Eqs  (12)  and  (17).  If 
the  flow  starts  with  frozen  Initial  con¬ 
ditions  It  is  known  to  remain  in  the  fro¬ 
zen  sta-e.  Hence,  the  solution  is  obtained 
from  the  simple  algebraic  equations  corres¬ 
ponding  to  the  frozen  flow  case.  If  rhe 
flow  starts  with  an  initial  state  which  is 
in  nonequilibrium,  the  similar  governing 
equations  (Eqs  (12)  and  (17))  have  to  be 
solved  using  the  specified  initial  values 
of  i}>,  and  n  for  the  given  parameters  X2, 

D  ard  i j .  However,  in  almost  all  practical 
situations  of  nozzle  flows,  the  flow  starts 
in  the  nozzle  reservoir  with  equilibrium 
conditions  and  remains  in  near-equilibrium 
up  to  the  nozzle  throat.1*8  This  will  be¬ 
come  evident  when  the  variation  of  (Ns) , 
along  the  nozzle  axis  is  discussed  in  the 
following  section.  Therefore,  in  the 
present  analy&io  the  flow  is  assumed  to  be 
in  an  equil-*hriura  condition  up  to  the 
throat.  In  this  case  the  function  41  “  'J* 


p  =  (const) (T) 5/?  =  (const)T 


(19) 


it 


*nd  the  Initial  value  of  t|/  cun  be  obtained 
from  the  equilibrium  oolution  {Eq  <21>)  for 
Value*  of  h  and  S0.  In  this  rather  in¬ 
direct  way  the  entropy  S0  also  appears  as  a 
barassater  in  the  nonequilibriua  solutions. 
It  is  noted  that  this  behavior  is  very  sim¬ 
ilar  to  that  of  the  dlssociatiocal  non- 
aquillbrium  case.8 

For  a  given  gao  (D  constant)  and  a  fam¬ 
ily  of  nozzle  shapes  (ij  constant)  the 
parameters  X2  and  SQ  nave  to  be  specified 
for  the  nonequilibrium  solutions .  However, 
this  two  parametric  dependence  car»  be  re¬ 
duced  to  a  single  one  by  Che  following 
transformation.  A  new  variable  £  is  de¬ 
fined  as  £  »  (S0  -  n) .  Then  Eq  (21) ,  which 
gives  the  equilibrium  solutions,  reduces  to 

5  -  -f-  -  loge  [^(l-e-*)]  ♦  f,  ♦  Sr 
6  (22) 


reservoir  conditions.  It  was  observed 
from  a  number  of  nozzle  computations  that 
the  velocity,  when  nondimensionalized  with 
a  velocity  of  u'0  “  (p„/p£)i,  does  not 
change  very  much  for  different  reservoir 
conditions.  Therefore,  the  expression  for 
(Na)3  is  rewritten  as 


<ss>, 


where  u 


Sri  n. 

is 


-Wi, 


-  (u'/u^). 


r  ' 

(1+1/iJ) 

The  constant  term 


(25) 


Ug/ud) (l+X/ij)  ia  ^9  (25)  can  be  included 
in  the  expression  for  X2  in  Eq  (17)  and 
the  remaining  terms  can  be  taken  as  the 
nonsimilar  function  (Ng).,  namely: 


Thus,  the  equilibrium  solutions  may  be  rep¬ 
resented  by  a  single  universal  curve  show¬ 
ing  the  variation  of  $  with  the  independent 
variable  £.  The  governing  similar  equa¬ 
tions  (Eq  (12)  and  (17))  for  the  nonequi¬ 
librium  case  reduce  to 

(5/2)d(log  fc)*d£+<j»[e*/(e4,-l)2]d<i>(l-f  )  «  0 
e  (23) 


d$  _  1 

d£  =  («s)2 


eX2  e[(i-i/ijU-W 


1/3 


-«] 


(24) 


where 


<«,>.  ■ 


(iPo) 


X/3L.  \  d-J./  / J 

'  U1 


(26) 


and 


x  »  x  +  ioge(ud/u^)(1+1/i’5) 

1  z  e  d  o  (27J 

Then  the  parameter  X2  defined  in  Eq  (26) 
becomes 

xi =  "  n-i/ij)s0]  ^28j 

By  letting  i  *  1.0  and  j  •>  2.0  (correspond¬ 
ing  to  the  hyperbolic  nozzle  case)  Eq  (26) 
is  rewritten  as: 


(N 


3/2 


(29) 


x2  =  [x2  -  (i-i/i j)soi 

The  nonequilibriun  solutions  depend  on 
a  single  parameter  x2  for  a  given  gas  and 
for  a  family  of  nozzle  shapes.  The  initial 
values  for  Che  noncqulllkrium  solutions  can 
be  easily  obtained  from  Eq  (22)  for  equi¬ 
librium  starting  conditions.  The  parameter 
X2  contains  the  effects  of  the  vibrational 
relaxation  time  constant  tv,  the  reservoir 
conditions  and  the  nozzle  shape.  The  non- 
similar  function  (N8)2  has  to  be  taken  into 
account  Jr.  order  to  obtain  an  exact  non- 
equilibrium  solution  from  Eqs  (23)  and  (24). 
The  properties  of  (t»3)2  and  a  method  of  its 
correlation  are  discussed  in  the  following 
section . 

Nona 1ml la r  Function  Ns 

The  function  (Ns)2  given  in  Eq  (24)  is 
a  function  of  H,  u  and  A.  Hence,  it  will 
have  different  values  for  different  res¬ 
ervoir  conditions  and  it  also  varies  along 
the  nozzle  axis.  The  velocity  u  ir.  the 
expression  for  ( N s 1  has  been  nondimension- 
alizod  by  ud  which  Is  independent  of  the 


A  typical  variation  of  (N3)  2  vich  area 
ratio  is  shown  in  Fig  1.  It  is  noted  that 
moving  upstream  from  the  nozzle  throat, 

(Kg) j  rapidly  tends  to  zero  since  H  and  u2 
both  ■+  0.  At  the  geometric  throat  the  rate 
of  change  of  area  with  x  goes  to  zero  for 
all  nozzles  with  co  area  discontinuity  at 


Figure  1.  Typical  Variation  of 

Function  (N„),  with  Area 
Ratio  for  Nitrogen 
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the  throat.  Therefore,  (Ns) s  has  an  in¬ 
determinate  fora  since  H  also  becomes  unity. 
Hovever,  it  can  be  shown  that  It  tends  to  a 
definite  Halt  t t  the  throat.  To  obtain 
this  Unit,  the  complete  expression  for 
(dn/dx)  given  below  must  be  considered. 


=.  1  M 


(=&) 


If  an  area  discontinuity  exists  at  the 
throat  (for  example  at  the  Juncture  of  two 
conical  nozzles)  the  rate  of  change  of  area 
has  a  finite  value  and  (dn/dx)  tends  to  in¬ 
finity  at  the  throat.  However,  nozzles 
used  in  practice  generally  have  no  area  dis¬ 
continuity  at  the  throat  so  that  the  func¬ 
tion  (Na)t  is  expected  to  behave  as  shown 
in  Fig  1.  The  function  (Ng) j  remains  pos¬ 
itive  even  for  M  <  1.0  in  the  upstream 
portion  of  the  nozzle  since  dA/dx  is  neg¬ 
ative  and  hence  the  negative  sign  iu  Eq 
(29)  applies  to  the  upstream  portion  from 
the  nozzle  throat. 

Vibrational  nonequilibrium  nozzle  flow 
quantities, starting  at  the  nozzle  throat, 
were  computed  for  nitrogen  for  a  number  of 
reservoir  conditions  by  using  the  computer 
program  of  Ref  11.  The  £  values  were  also 
computed  using  the  equation  5  “  (S0-n)  = 

(Sc  +  logep) .  The  reservoir  entropy  was 
computed  from  Eq  (21)  where  the  reference 
entropy  Sr  was  taken  equal  to  15.2  for 
nitrogen  which  makes  the  entropy  values  the 
sane  as  those  computed  in  Ref  12.  The  (Hs) j 
values  were  also  computed  for  several  res¬ 
ervoir  conditions,  using  Eq  (26),  and  are 
plotted  in  Fig  2  with  (5^/5)  as  the  vari¬ 
able.  Although  all  the  curves  correlate 
fairly  well  immediately  downstream  of  the 
nozzle  throat,  a  significant  temperature 
effect  shows  up  for  §*/£  values  greater 
than  about  1.2.  In  an  effort  to  obtain  a 
better  correlation  the  following  approach 
was  taken.  In  the  mass  flow  correlation 
analysis  (Appendix  B)  it  is  noticed  that 
the  nondlnensionol  mass  flow  (p^i/o^u^)  Is 
also  slightly  temperature  dependent.  There¬ 
fore,  the  (Ns) t  values  were  multiplied  by 


che  corresponding  mase  flow  ratio  as  given 
below 

«S  =  (NsMe*u>oUo36'°  (31) 

where  the  exponent  6.0  vos  determined  by 
equating  the  ha  values  corresponding  to  the 
maximum  and  minimum  valuea  of  (NB),  given 
in  Fig  2.  The  Na  values  computed  "from 
Eq  (31)  are  shown  in  Fig  3  for  the  same 
reservoir  conditions  used  in  Fig  2.  All 
the  values  for  different  reservoir  con¬ 
ditions  correlate  very  well  and  can  be 
represented  by  a  mean  curve  as  shown  by  the 
dotted  line  la  Fig  3.  The  maximum  dis¬ 
crepancy  of  the  actual  values  from  the  mean 
curve  is  within  a  few  percent.  Tbe  mean 
curve  can  be  also  represented  by  a  simple 
analytical  equation  of  the  form: 


Ns  =  0.37  -  0.32  (2.0  - 


U*/£)  <  2.0 


Ns  *  0.37  for  tejtf  >  2.0 

The  Ns  valuea  were  also  computed  with 
different  hyperbolic  nozzle  shapes  (L  “ 

0.5  to  2.0)  as  well  as  combined  hyper¬ 
bolic  (L  “  0.5  to  2.0)  and  conical  (L  ■*  1 
to  3)  shapes.  It  waB  found  that  the  dif¬ 
ferences  iu  the  Ns  values  were  within  the 
accuracy  of  the  correlation  shown  in  Fig  3. 
The  factor  [Piui,*pouol6*9  uas  also  in¬ 
cluded  in  the  expression  so  that  its 
effect  is  properly  taken  into  account. 

After  including  che  mass  flow  factor,  Che 
final  expressions  for  Ng ,  1  and  x  ate. 
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Figure  2.  Variation  of  Function 

(Ha),  with  the  Parameter 
(5*/5)  for  Nitrogen 


Figure  3.  Correlation  of  Function 
Na  with  the  Parameter 

/  r  /  r\  ui.. _ 


The  two  similar  governing  equations 
(Eqs  (23)  and  (24)  with  Ns  given  by  Eq  (32) 
were  solved  by  a  fourth-order  Runge-Kutta 
technique  with  ij  *•  2.0  and  D  -  14.7  which 
corresponds  to  nitrogen.'3  The  starting 
values  of  (i  and  $  were  obtained  for  a  given 
£  from  the  equilibrium  solution  ($  “  $) 
given  by  Eq  (22) .  The  reference  entropy  Sr 
in  Eq  (22)  was  cakcn  equal  to  15.2.  The 
factors  f j ,  f  ,  ,  and  f,  are  all  assumed  to 
be  zero  (see  Appendix  A).  The  starting 
values  of  £  for  different  x  values  were 
selected  in  such  a  way  that  the  solution 
always  starts  with  equilibrium  conditions. 
The  equilibrium  solution  obtained  from 
Sq  (22)  is  shown  in  rig  4  and  i3  repre¬ 
sented  by  a  single  universal  curve  since  it 
is  independent  of  Xt  ij  and  D.  Also,  in 
Fig  4  a  typical  similar  nonequilibrlum 
solution  is  compared  with  an  exact  solution 
obtained  from  the  computer  program  of  Ref 
11.  The  comparison  is  considered  ro  be 
very  good. 

A  series  of  solutions  for  different 
values  of  x  for  nitrogen  (D  »  14.7)  are 
shown  in  Fig  S.  The  vibrational  tempera¬ 
ture  function  -J  is  seen  to  follow  the  trans¬ 
lational  temperature  function  £  very  closely 
for  awhile,  the  extent  of  which  depends  on 
X.  and  then  diverges  rather  suddenly  and 
reaches  a  constant  value:  this  corresponds 
to  the  freezing  of  the  vibrational  energy 
node.  The  translational  temperature 


( 


Figure  4.  '’onparison  of  Present  Similar 
Solutions  with  an  Exact 
Solution  <x  *  4.4,  D  “  14.7, 
ij  -  2.0 ) 


Figure  5.  Similar  Solutions  for 

Vibrational  Nonequilibrium 
Nitrogen  Flow  (D  *>  14.7) 


function  i[»  increases  monotonically  as  £ 
decreases.  The  equilibrium  solution  shown 
in  Fig  5  is  also,  given  by  the  envelope  of 
all  the  nonequilibrium  solutions.  The  con¬ 
stant  D  “  14.7  for  nitrogen  was  obtained 
from  the  correlation  of  experimental  values 
over  a  certain  temperature  range.9  How¬ 
ever,  1)  can  have  slightly  different  values 
for  the  same  gas  over  different  ranges  of 
tenperatures.  Therefore,  similar  solutions 
were  also  computed  with  D  °  13.5  and  15.5 
and  are  shown  in  Figs  6  and  7  respectively. 
It  would  be  an  easy  matter  to  interpolate 
between  these  solutions  for  slightly  dif¬ 
ferent  values  of  D.  The  frozen  vibrational 
temperature  function  $£  depends  on  x  only 
and  the  variation  of  with  x  for  all  the 
three  values  of  D,  is  shown  in  Fig  8. 


«  6  8  KJ  12  W  16  18  CO  22 


Figure  6.  Similar  Solutions  for 

Vibrational  Nonequilibrlum 
Nitrogen  Flow  (D  m  13.5) 

V.  THE  PARAMETERS  g  AND  x 

It  is  shown  in  this  analysis  that  the 
nonequilibrium  s<milar  solutions  depend  on 
two  general  parameters  £  and  x*  In  order 
to  use  the  similar  solutions  presented  in 
this  report  the  parameters  ?  and  x  should 
be  known  in  terms  of  the  initial  and  bound¬ 
ary  values.  Therefore,  the  functional  de¬ 
pendence  of  £  and  x  are  considered  in  this 
section . 
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Figure  7.  Similar  Soli. Cions  for 

Vibrational  Nonequilibrium 
Nitrogen  Flow  (D  »  15.5) 


Figure  8.  Frozen  Vibrational  Temperature 
Function  with  Similarity 
Parameter  x  for  Nitrogen 


Figure  10.  Velocity  Function 

Correlation  with  Area 
Ratio  for  Nitrogen 


[0.5  -  0.31(1  +  logi0A)-2]  (37) 


Faramecer  S 

The  parameter  C  is  defined  as 


This  parameter  is  not  only  a  function  of  the 
reservoir  and  nozzle  throat  conditions  but 
also  a  function  of  velocity,  which  is  un¬ 
known.  Therefore,  the  velocity  ratio  u'/uj 
was  computed  for  several  sets  of  reservoir 
conditions  and  its  variation  with  nozzle 
area  ratio  is  shown  in  Fig  9.  It  is  ob¬ 
served  that  this  type  of  correlation  still 
results  in  a  significant  amount  of  reservoir 
temperature  effect  which  is  similar  to  that 
noticed  in  the  correlation  of  (Ns)j  shown 
in  Fig  2.  Therefore,  eacn  velocity  ratio 
was  multiplied  by  the  corresponding  mass 
flow  ratio  and  replotted  In  Fig  10.  All  the 
computed  values  now  correlcte  very  well  and 
a  mean  curve  can  be  orawn  through  the  points 
as  shown.  This  mean  curve  can  be  repre¬ 
sented  by  a  simple  expression  of  the  type 


The  exponent  4.0  in  Eq  (37)  was  determined 
by  using  considerations  similar  to  those 
used  previously  in  the  correlation  of  Ns. 
By  substituting  for  Sc,  u'/u*  and  p*  from 
Eqs  (21),  (37)  and  (B2) ,  respectively,  an 
expression  for  £,  after  some  algebraic 
manipulations,  may  be  given  as 


log^ 


(l-e"Vo) 


X 


A[0.5  -  0.31(1+1ogw 


+  S_ 


(38) 


It  Is  noted  that  £  depends  on  only  two 
parameters,  namely,  the  ares  ratio  A  and 
the  reservoir  temperat-ire  function  0o.  The 
parameter  £  does  not  depend  on  either  the 
reference  density  pj  or  the  reservoir  den¬ 
sity  Pg.  The  parameter  £  is  now  expressed 
in  terms  of  the  initial  and  boundary  values 
only  and  hence  can  be  readily  computed 
for  any  given  conditions. 
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Parameter  x 


•1/iJ,.,  {  1/2+1/ IJ  ) .  (1-1/ij) 

"o  *o 


(41) 


where  the  uass  flow  correlation  expression 
given  in  Appendix  B  has  been  used  to  elira- 
inate  the  mass  flow  tern  in  X.  With  Eqs 
(40)  and  (41) ,  the  parameter  X  reduces  to 


(6+1/1 


(i»-2.5/ij) 


-d-l/i 


(42) 


It  is  noted  that  the  terns  (R9v/m)l,  L  end 
C/po  have  the  dimensions  of  velocity,  length 
and  tine,  respectively.  Vo  an<*  Sr  are  non" 
dinensional  nunbers.  If  pj,  is  ueasured  in 
atnospheres  then  the  C  units  will  be  ata- 
s«c.  The  paraneter  x  X-8  independent  of  A. 
For  i  given  gas  x  depends  on  only  p^,  L  and 
Vo  »ince  0V,  C  and  5r  arc  all  constan.  . 

The  functions  $  at d  V  as  well  ns  the  gen¬ 
eral  paraneters  £  and  X  do  net  depend  on 
any  of  the  reference  values  that  are  used 
for  nondioensiunallzlng  the  governing  equa¬ 
tions.  Hence,  the  reference  values  can  be 
chosen  arbitrarily. 


In  Kef  6  the  paraneter  p£L  was  used  as 
a  correlating  paraneter.  Also,  in  Ref  7  a 
paraneter  (p^L  y^i/C)  was  deduced  by  non- 
dltaensionali ring  the  sudden  freezing  cri¬ 
teria  with  reservoir  values.  These  two 
parameters  can  be  deduced  from  the  present 
general  paraaeter  *.  To  obtain  the  param¬ 
eter  PqL  the  entire  effect  of  Vo  ^-a8  to  be 
neglected.  The  other  parameter  can  be  ob¬ 
tained  from  x  hy  neglecting  a  portion  of 
the  effect  of  v0*  This  neglect  of  the  V0 


effect  Is  the  reason  why  the  frozen  vi¬ 
brational  temperature  presented  in  these 
analyses6)7  depends  on  the  reservoir  tem¬ 
perature  in  addition  to  the  approximate 
paraneters.  Furthermore,  these  analyses 
predict  only  the  approximate  frozen  vi¬ 
brational  temperature  and  do  not  provide 
the  flow  quantities  in  the  nonequilibrium 
region. 

Range  of  the  Applicability 
of  the  Paraneter  x _ 

The  general  correlating  parameter  X»  as 
noted  before,  depends  on  Pq,  I,  and  Vo  ^or  8 
given  gas  .  The  variation  of  with  Vo 
for  a  constant  x  can  be  computed  from 
Eq  (42)  and  is  shown  in  Fig  11  for  a  number 
of  x  values  with  ij  ■=  2.0.  The  variation 
of  Pq  (L  <*  1.0)  with  Vo  f°r  3  constant 
equilibrium  mole  fraction  of  O.i  is  also 
shown;  this  curve  represents  approximately 
the  high  temperature  limit  beyond  which  dls 
dissociation  becomes  appreciable  and  the 
dissociation  relaxation  nay  have  to  be  con¬ 
sidered  in  addition  to  the  vibrational  re¬ 
laxation  phenomenon.  It  can  be  inferred 
from  Fig  il  that  all  of  the  practically 
feasible  reservoir  conditions  (a  maximum 
pressure  of  1000  atmospheres  and  a  tempera¬ 
ture  range  frou  2000#K  to  8000°K)  are  cov¬ 
ered  by  x  values  between  1  to  8,  which  is 
rather  a  narrew  range  •  jmpared  to  the  range 
of  x  values  presented  in  Fig  5.  It  is  also 
observed  in  Fig  5  that  the  nonequilibrium 
solution ,  for  a  given  x>  departs  from  the 
equilibrium  solution  at  a  certain  maximum 
Value  of  V;  these  maximum  V  values  are  also 
plotted  in  Fig  11.  For  reservoir  condi¬ 
tions,  which  fall  above  this  line,  the 
solutions  start  with  equilibrium  conditions 
and  can  be  obtained  from  the  present  sim¬ 
ilar  solutions.  For  reservoir  conditions 
which  fall  well  below  this  line  the  flows 
can  be  taken  as  completely  frozen  in  the 
entire  nozzle.  In  a  narrow  region  just 
below  the  equilibrium  limit  line  the  flow 
will  be  in  the  nonequilibrium  state  and  the 


■d? 

£ 


Figure  11.  Range  of  Applicability 

of  the  Parameter  x  (ij  "  2.0, 
C  »  1.715  x  10”11  atm-sec., 

Po  in  atm.,  L  in  cm) 
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Solutions  have  to  be  obtained  by  starting 
with  reservoir  conditions  as  the  initial 
values.  Furthermore,  the  function  Ns  up¬ 
stream  of  the  throat  has  to  be  also  in¬ 
cluded  in  the  solutions. 

Effect  of  Starting  Value 
£#  on  Similar  Solutions 


It  is  observed  in  Fig  5  that,  for  a 
given  value  of  X»  there  is  a  range  of 
(hence  £*)  values  that  could  be  used  as  the 
initial  starting  values.  The  upper  limit 
is  given  by  the  equilibrium  limit  line  in 
Pig  11  and  the  lower  limit  is  given  by  the 
dissociation  limit  as  shown  in  Fig  11.  For 
example,  the  range  of  values  is  approx¬ 
imately  1.1  to  0.6  for  X  “  8.0  and  C.65  to 
0.45  for  x  “  4.0.  Hence,  for  a  given  x« 
the  nonequilibrium  solutions  could  be 
started  with  any  of  a  range  of  5*  values. 
This  would  result  in  slightly  different 
values  of  Ns  and,  therefore,  different  non- 
equilibrium  solutions  for  the  same  value  of 
X>  To  examine  this  point  more  closely, 
nonequilibrium  solutions  were  obtained  for 
different  £A  values  but  with  the  same  x 
value  and  it  was  found  that,  for  the  range 
of  allowable  starting  values  mentioned  pre¬ 
viously,  the  differences  in  the  nonaqui- 
librium  solutions  varied  only  a  few  percent 
and  this  error  is  within  the  accuracy  of 
the  correlation  of  Na  shown  in  Fig  3. 


VI.  CONCLUSIONS 


Eased  on  the  present  analysis  the  fol¬ 
lowing  conclusions  are  reached: 

1.  Similar  solutions  for  vibrational  non¬ 
equilibrium  nozzle  flow  problems  can  be 
obtained  by  using  the  new  similarity  param¬ 
eter  £. 

2.  The  similar  solutions  presented  can  be 
used  over  a  wide  range  of  practicable  com¬ 
binations  of  initial  conditions  and  nozzle 
scale  parameters. 

J.  The  vibrational  equ* librium  solutions 
depend  on  the  one  parameter  €  only  and  the 
nonequilibrium  solutions  depend  on  two 
parameters  £  and  x*  The  frozen  vibrational 
temperature  depends  on  X  only. 


vibratlonally  excited  through  an  infinite 
number  of  equally  spaced  energy  levels.  In 
the  real  situation,  the  excited  molecules 
dissociate  when  the  vibration  energy  level 
corresponding  to  the  dissociation  energy  la 
reached.  Therefore,  in  the  cut-off  har¬ 
monic  oscillator  approximation  the  vibra¬ 
tion  energy  is  considered  only  up  to  the 
dissociation  limit.  Since  the  energy  lev¬ 
els  are  assumed  to  be  equally  spaced,  the 
number  H  of  allowable  energy  levels,  when 
applying  this  approximation  is  given  by 

N  ■  W 

The  factors  whlcn  take  into  account  the 
effect  of  the  cut-off  harmonic  oscillator 
approximation  are  given  belo:;: 


nV*V-i)2 


(Al) 


_(e*-e*)(eH*-l}(eN*-l)  J  (A2) 


f 

3 


|joge(1-e"K*) 


(A3) 


It  can  be  easily  shown  that  these  factors 
tend  to  zero  as  N  This  limit  corre¬ 

sponds  to  the  simple  harmonic  oscillator 
approximation.  Since  no  other  expression 
in  the  governing  equations  contains  the 
parameter  N,  the  governing  equations  for  a 
simple  harmonic  oscillator  approximation 
can  be  obtained  by  simply  assuming  these 
factors  are  equal  to  zero.  Furthermore,  it 
can  be  shown  that  for  nitrogen,  with  N  -  34, 
the  correction  factors  are  all  very  small 
even  for  temperatures  as  high  as  7000*K. 

For  example,  at  Tj  “  T^  -  6000*K  and  N  »  34, 

f,  -  2.28  x  10_t,  f2  -  0,  f3  »  -8.7  x  10~*. 

Thus,  for  the  problem  under  consideration, 
the  simple  harmonic  oscillator  model  is 
more  than  adequate. 


APPENDIX  B 


4.  The  parameters  £  and  X  serve  as  univer¬ 
sal  correlating  parameters  since  they  con¬ 
tain  all  the  parameters  of  the  problem. 

5.  The  present  similarity  transformation 
of  the  governing  equations  affords  a  better 
insight  of  the  parametric  dependence  in 
this  problem  and  should  be  of  interest  to 
tb -oruticiens  as  well  as  experimentalists. 


APPENDIX  A 

The  Cut-Off  Harmonic 
Oscillator  Approximation 

The  simple  harmonic  oscillator  model 
assumes  the  diatomic  molecule  may  be 


Correlation  of  Kass 
Plow  and  Throat  Density 

The  general  correlating  parameters  £ 
and  x  contain  the  mass  flow  p*u^  and 
critical  throat  density  p*  which  have  to 
be  obtained  in  a  separate  computation.  If 
the  flow  at  the  throat  is  in  nonequillbrium, 
then  the  quantities  p*o^  and  p^  have  to  be 
computed  by  a  trial  and  error  procedure. 
However,  in  the  present  analysis  the  flow 
is  considered  to  be  in  vibrational  equi¬ 
librium  up  to  the  throat  for  which  situation 
Che  mass  flow  and  the  throat  density  can  be 
obtained  by  a  set  of  algebraic  equations.  In 
the  present  case,  the  mass  flow  values  were 
computed  for  nitrogen  over  a  wide  range  of 
reservoir  conditions  by  using  the  computer 


program  of  Ref  12  and  ore  presented  in  Fig  2.  Anderson,  J.D.,  "A  Tine-Dependent  Annl- 
Bl,  The  nondimensional  mass  flow  (pAu|)/  ysie  for  Vibrational  and  Chemical  Nonequi- 
(Pouo^  is  independent  of  reservoir  presoure  librium  Nozzle  Flows",  AIAA  Journal,  Vol  8, 
but  slightly  dependent  on  reservoir  tern-  No.  3,  Mar  1970,  pp.  545-550. 
perature.  It  can  be  represented  by  a 

linear  equation  of  the  type  3.  Sto.llery,  J.L.  &  Smith,  J.E.,  "A  Note 

on  the  /ariation  of  Vibrational  Temperature 
Along  a  Nozzle,  Journal  Fluid  Mcch.,  Vol  13 

pX  -  MpW  (B1)  1962‘  pp-  225-235- 

4.  Stollery,  J.L.  &  Park,  C.  ,  "Couputer 

wbere  Solutions  to  the  Problem  of  Vibrat1* onal 

g  ,  Relaxation  in  Hypersonic  Nozzle  Floys", 

k,  «  [0.689  -  6.3  x  10“  T  (°K)3  Journal  Fluid  Mach. .  Vol  l.\  1964,  pp.  113- 

123. 


A  similar  correlation  for  the  nondimension- 
al  throat  density  p^/p'  is  also  shown  in 
Fig  Bl.  It  i«  alR«  independent  of  reser¬ 
voir  pressure  and  can  oe  represented  by  a 
linear  equation  of  the  type 

p*  ■  Mo  <32> 

where 

k,  -  [0.634  -  2.33  x  10"6  r  (°K)3 


Figure  Bl.  Correlation  of  Throa 
Density  and  Mass  Flo* 
Ratios  for  Nitrogen 
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